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Abstract. For multiquadric interpolation it's well known that the most powerful error bound 
is the so-called exponential-type error bound. It's of the form \f(x) - s{x)\ < CX^ d \\f\\ h where 

< A < 1 is a constant, C is a constant, d is the fill distance which roughly speaking measures the 
spacing of the data points, s(x) is the interpolating function of f(x), and \\f\\h denotes the norm off 
induced by the multiquadric or inverse multiquadric. This error bound for x £ R n was put forward 
by Madych and Nelson in 1992 and converges to zero very fast as d — > 0. The drawback is that 
both C and A get worse as n — ► oo. In particular, A — ► l~ very fast as n — > oo. In this paper, we 
raise an error bound of the form \f(x)-s(x)\ < C'Vd\ n / d \\f\\ h where A' can be independent of the 
dimension n and < A' < A < 1. Moreover, C is only slightly different from C. What's noteworthy 
is that both C and A' can be computed without slight difficulty. 
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1 Introduction 

We begin with some fundamental knowledge about radial basis functions. Let h be a continu- 
ous function in R n and be strictly conditionally positive definite of order to. Given data points 
(xj,fj),j = 1,...,N, where X = {x\, . . . ,xjv} is a subset of points in R n and the /js are real or 
complex numbers, the so-called h spline interpolant of these data is the function s defined by 

N 

s(x) =v{x) +^2cjh(x- Xj), (1) 

where p(x) is a polynomial in P^ l _ 1 , the space of n-variable polynomials of degree not more than 
to — 1, and the c/s are chosen so that 

N 

Y, c il(. x j) = ° ( 2 ) 
i=i 
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for all polynomials q in P^_i and 

N 

P( x i) + ^c j h(x l - Xj) = f i} i=l,...,N. (3) 

3=1 

It is well known that the linear system (2) and (3) has a unique solution if X is a determining set for 
Pm-i an d ^ ls strictly conditionally positive definite. Further details can be found in [7]. Therefore 
in our case the interpolant s(x) is well defined. 

Here, X is a determining set for Pj^-i means that X does not lie in the zero set of any nontrivial 
member of P^_ x - 

In this paper h is defined by 

h(x) :=r(-|)(c 2 + |af)§, /3€i?\27V> , c> 0. (4) 

where \x\ is the Euclidean norm of x, T is the classical gamma function and c are constants. 
The function h is called multiquadric or inverse multiquadric, respectively, depending on (3 > 0, or 
(3 < 0. 

In [3] Madych and Nelson raise the famous exponential-type error bound for scattered data inter- 
polation of multiquadric and inverse multiquadric, as mentioned in the abstract. The computation 
of C and the crucial constant A is solved by Luh in [5j . This brings the error bound forward to the 
possibility of application. However, in practical application, too many data points may result in a 
big condition number for solving the linear system (2) and (3). A good solution is to improve the 
error bound so that a satisfactory error estimate can be reached before ill-conditioning occurs. This 
is the purpose of this paper. 



1.1 Some Fundamental Knowledge about Polynomials 

Let P™ be the space of polynomials of degree not exceeding I in n variables and E C R n be compact. 

It's well known that dimPJ 1 = ( n ^ ' ^. In this section we will denote dimPJ 1 by N. It's also 

well known that if x\, . . . , arjy G E do not lie on the zero set of any nontrivial q € P, n , there exist 
Lagrange polynomials U of degree / defined by h(xj) — Sij, 1 < i, j < N . This guarantees that 
for any / e C(E), (Uif)(x) := YliLi f{ x i)k{x) is its interpolating polynomial. It's easily seen that 
Hj(p) = p for all p £ Pp and hence the mapping n; : C(E) — ► P" is a projection. Let 

AT 

j|H ; || i-max^l^)!. 

i=l 

It's easy to show that for any p £ P™, 

Iblloo < ||n ; || max . \p(xi)\ 

l<t<N 

if the domain of p is E. In this paper we will discuss mainly the case E — T n , an n-dimensional 
simplex whose definition can be found in [B]. 

Now we are going to adopt barycentric coordinates. Let u,, l<i<n + l, be the vertices of 
T n . Then any x £ T n can be written as a convex combination of the vertices: 

n+1 n+1 

X = ^ C i V ii X! C,: = 1 ' Ci ~ - 

i=l i=l 



2 



We define "equally spaced" points of degree / to be those points whose barycentric coordinates 
are of the form 



ra+1 



(ki/l, k%/l, . . . , k n+ i/l), ki nonnegative integers with fc, = I, 

i=l 

It's obvious that the number of such points is exactly N = dimP™ . Moreover, by pQ, we know that 
equally spaced points form a determining set for P". We will need the following lemma which we 
cite from pQ. 

Lemma 1.1 For the above equally spaced points ||II/| < ( ^ . ^ ]. Moreover, as n — > oo, ||II;|| — > 



21 - 1 
i 



Now, we are going to prove a lemma which plays a crucial role in our construction of the error 
bound. 

Lemma 1.2 Let Q C R n be a simplex and Y be the set of equally spaced points of degree I in Q. 
Then, for any point x in Q, there is a measure a supported on Y such that 

p{y)do-(y) =p{x) 

for all p in P™, and 

21-1 



d\o-\(y) < 

Proof. Let Y — {yi, . . . ,Vn} be the set of equally spaced points of degree I in Q. Denote P™ 
by V. For any x £ Q, let 8 X be the point-evaluation functional. Define T : V — > T(V) C R N 
by T(v) = (S yz (v)) yzeY . Then T is injective. Define $ on T(V) by ip(w) = S^T^w). By the 
Hahn-Banach theorem, ip has a norm-preserving extension ipext to P w . By the Riesz representation 
theorem, each linear functional on R N can be represented by the inner product with a fixed vector. 
Thus, there exists z £ R N with 

N 

4>ext(w) = ^2,Z j Wj 
3=1 

and ||x|bjjjv)* = HV'extll- If we adopt the Zoo-norm on R N , the dual norm will be the ^-norm. Thus 

ll^ll^), =\\z\\i = 11^**11 = M = \\S*T- 1 \\. 

Now, for any p £ V, by setting w — T(jp), we have 

N N 

8 x (p) = 5 x (T~ 1 w) = tf)(w) = 4>ext{w) = ^ZjWj = ZjS y .(p). 

3=1 3=1 

This gives 

N 

P( X ) =^2 Z 3P(V3) ( 5 ) 
3=1 

where \zA + ■■■ + \z N \ = II^T- 1 !!. 
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Note that 



Therefore 



\6 X T II = sup 



< 



w e T(F) HI*" 
IM 



sup 



wET{V) ITOIU* 



< sup W*>l 



< sup 



pey max^i,...,^ \p{Vj)\ 

III || max 3= i ; ... jA r |p(2/j)l 



pey max j= i,...,iv |p(l/j)| 
Hll 

2; - 1 



21 - 1 
Z 



|*l| + --- + M < 

and our lemma follows immediately by (5) and letting c{{yj}) = Zj, j = 1, . . . , N. 



1.2 Some Fundamental Knowledge for Interpolation 

In our development of the error estimate some basic ingredients should be introduced. We emphasize 
on two objects:first, the function space where the interpolating and interpolated functions come 
from;second, an important lemma introduced by Gelfand and Vilenkin in [2], but modified by Madych 
and Nelson in [8 , which reflects deeply the features of a radial basis function. 

Let V denote the space of complex-valued functions on R n that are compactly supported and 
infinitely diffcrcntiablc. The Fourier transform of a function <f> in D is 

4>(0 = J e- i<x ' s> 4>(x)dx. 

A continuous function h is called conditionally positive definite of order m if 

J h(x)(f>(x) * 4>{x)dx > 

holds whenever <f> — p(D)ip with ip in T> and p(D) a linear homogeneous constant coefficient differ- 
ential operator of order m. Here <f> = <j)(—x) and * denotes the convolution product 

<f>x * (f>i(t) = J (j>i(x)<f>2(t — x)dx. 

This definition, introduced in [S], looks complicated. However it's equivalent to the generally used 
definition. For further details we refer reader to [7] and [5|. 
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Now we have come to the afore-mentioned crucial lemma. It's as follows. If h is a continuous 
conditionally positive definite function of order m, the Fourier transform of h uniquely determines 
a positive Borel measure /i on i?™\{0} and constants a 7 , 17I = 2m as follows: For all tp € V 



1 



h(x)^(x)dx = I E Dr $(o)^\M€) ( 6 ) 

h\<2m ' J 

E ^(o)^. 

|7|<2m 

where for every choice of complex numbers c a , \ a\ = m, 

E E a "+/ 3c "^' - °- 

|a|— m |/3|— m 

Here x is a function in T> such that 1 — x(£) has a zero of order 2m + 1 at £ = 0; both of the integrals 
/o<|£|<i l£| 2m ^/' t (£)) Jjfi>i ^A*(0 are nn he. The choice of \ affects the value of the coefficients a 7 for 
I7I < 2m. 

Our error estimate is based on a function space called native space whose characterizations can 
be found in [3], [4] and [11]. Although there are simpler expressions for the function space, we still 
adopt the one used by Madych and Nelson in [8] to show the author's respect for them. First, let 

V m = j<£ e V : J x a (/)(x)dx = for all \a\ < m| , 

Then the native space Ch, m is the class of those continuous functions / which satisfy 



f(x)4>(x)dx 



<c(/) / h(x-y)4>{x)<l>(y)dxdy} (7) 



1/2 



for some constant c(/) and all (j> in V m . If / G Ch, m , let ||/|U denote the smallest constant c(/) for 
which (7) is true. Here is a semi- norm and Ch. m is a semi-Hilbert space; in the case m = it 
is a norm and a Hilbert space respectively. 



2 Main Results 

It's well known that the function h defined in (4) is conditionally positive definite of order m = if 
(3 < 0, and m = \^~\ if /3 > 0. With this in mind we have the following lemma. 

Lemma 2.1 Let h be as in (4-) and m be its order of conditional positive definiteness. There exists 
a positive constant p such that 

f \£\ l dn(0 < (V2)" W1 • • na n ■ cP~ l ■ A • P l ■ 11 (8) 

for all integer I > 2m + 2 where fi is defined in (6), a n denotes the volume of the unit ball in R n , c 
is as in (4), and Aq is a positive constant. 
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Proof. Let K v denote a modified Bessel function of the second kind. Then 

n + P 



2c 

+13 



L 
- L 

= 2 - n/2 (£f 2 Jjzt^^mw 

Tt+g+i n+i o_. f 00 r , , n — L 

= 2 2 7r 2 na n c p / — dr where i = I H 

7o e r 2 



Note that if /3 < 0, then m = and Z > 2m + 2 = 2. This implies Z' > 0. If /? > 0, then m = |~§] 
and Z > 2m + 2 = 2[|] + 2. This implies Z' > 0. In any case Z' > 0. 

Our proof then proceeds in three cases. Let I" = which is the smallest integer greater than 
or equal to I'. 

Casel. Assume I" > I. Let Z" = I + s. Then 



/■°o J r°° r l" 

/ — dr< / — dr = l"\ = (Z + s)(Z + s-l) •••(« + 1)Z! 
Jo e Jo e 

and 



/ ~^-dr< / -_rf r = (Z" + l)! = (Z + s + l)(Z + s)---(Z + 2)(Z + l)Z!. 
Jo e Jo e 



Note that 

(Z + s + l)(Z + s)---(/ + 2) _ Z + s + 1 
{l + s)(l + s- 1) ■■■(/ + 1) ~ Z + l 

(i)Assume /3 < 0. Then m = and Z > 2. This gives 

Z + s + 1 3 + s 
Z + l ~ 3 

Let p= 2±». Then 



oo r+i 

^rfr <A 0P ' +1 (Z + 1)! 
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if J °° ^-dr < A p l ll. The smallest I" is l% = 2 + s. Now, 



r°° r l ° 

/ — dr = Z£! = (2 + s)(2 + s-l)---(3)f! where 1 = 2 
Jo eT 



(2 + S )(2 + r l)-(3) - ^ = 2) 



(2 + s)(2 + s- l)---(3) 



= Ao/5 2 2! where Ao = 

P 2 



It follows that J °° ^dr < A p J Z! for all I > 2. 

(ii) Assume > 0. Then m = [f ] and / > 2m + 2. This gives 



l + s + l 2m + 3 + s 

^ ~ " — 1 



l + l ~ 2m + 3 2m + 3 



r i" + l 



-dr < A p l+1 (l + iy. 



& 

if J °° ^p-dr < A p l l\. The smallest I" is Z£ = 2m + 2 + s when I = 2m + 2. Now, 
r°° r l 'd 
Jo e 

= IqI = (2m + 2 + s)(2m+ 1 + s) • • • (2m + 3) (2m + 2)! 
(2m + 2 + .s)(2m+l + S )...(2m + 3) 2m+2( „ 



. 2m +2^ „m , a (2m + 2 + s)(2m+ 1 + s) • • • (2m + 3) 
= A p 2m+2 (2m + 2)! where A = * A ^ 2m+2 ; i '- 



It follows that f °° ^dr < Aop'Z! for all I > 2m + 2. 

Case2. Assume Z" < L Let I" = 1 — s where s > 0. Then 

/•oo r ,00 /" 

/ — dr < / — dr = l"\ = (l- s)\ 
Jo e r ~ Jo e r 



l(l-l)---(l- s + l) 



and 



00 r ;'+i /-oo r ;"+i 

dr < I dr 

e r 



Note that 



= (r + l)! = (i -a + l)l = (i + 1)r . 1 (; _ a + 2) .( t+ l) ! . 



(i + 1)/ • • • (Z - s + 2) ' 1(1 - 1) • • • (I - s + 1) 
1(1 - 1) ■ ■ ■ (I - s + 1) 
(J + l)/---(i-« + 2) 
Z-s + 1 

m 
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(i) Assume (5 < 0, then m = and I > 2. Since Z" = I — s > 1 holds for all / > 2, it must be that 
s = 1. Thus 

Lzl + i = i _ _L_ = i _ JL < i /- or all l>2. 
l + l l + l l + l ~ J 

Let p = 1. Then 



1 



-dr< A 0( o' +1 (Z + l)! i/ / — dr<A p'Z!. 



/o e- 

1" ?o 7" 



The smallest Z" is Zg = Z - s = 2 - s. Now 



/ 

Jo 



e' 



dr = #! = (2 -«)! = !! = 1 



-Z! where I = 2 
2 

Aqp'Z! where A = — . 



It follows that L°° ^rdr < A p l ll for all Z > 2. 



o c 

(ii)Assume (3 > 0. Then m = |~§] and Z > 2m + 2. This gives 

s + 1 s 

= 1 - -, < 1. 



l+l l+l 

Let p=l. Then f °° 4^ dr ^ A o/°' +1 G + 1)! if / °° 4^~ dr ^ A o/^ ! - The smallest Z is Z = 2m + 2. 
Hence the smallest Z" is Z ' = Zo — s = 2m + 2 — s. Now, 

r°° Jn' 



r°° r l o 

/ — dr = Z£! = (2m + 2-s) = (Z -s)! 
Jo e 



(io)! 



Zo(Zo- 1) ■ ■ ■ (io - « + 1) 
= A 0( o'°Z ! w/iere A = 



(2m + 2) (2m + 1) • • • (2m - s + 3) ' 
It follows that J °° ^rdr < A p l l\ for all Z > 2m + 2. 
Case3. Assume I" = I. Then 

r°° J.' fOO I" roo i'+l 

/ _dr< / — dr = Z! and / dr < (I + 1)!. 



Jo 



Let p = 1. Then 

r°° r l' 

—dr < A p l l\ 
/o e 

for all Z where A = 1. The lemma follows now immediately from the three cases. 

Remark. In the preceding lemma the constant p plays an important role in our error estimate. 
Moreover, the constant A also means something in the error bound. Therefore we should make a 
thorough clarification for the two constants. It can be shown easily that 
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(a) I" > I if and only if n - > 3; 

(b) I" < I if and only if n - < 1; 

(c) I" = I if and only if 1 < n - < 3 

where I" and I are as in the proof of the lemma. 

The two constants can then be determined as follows. 

(a) n-0>3. Let s = [^f^] . Then 

(i) if < 0, then p = 2±£ and A = (2+')a+«)-3 . 

(ii) if > 0, then p = 1 + and A = (2™+2+0(2m+i+»)---(2m+3) where m = ^ . 

(b) n - < 1. Let s = - P^f^l . Then 

(i) if < 0, then p = 1 and A = ±; 

(ii) if /? > 0, then p = 1 and A = (2m+2)(2m+1) ... (2m _ s+3) i where m = ff ]. 

(c) 1< n - /3 < 3. Then p = 1 and A = 1. 

Before introducing our main result, we need the following lemma. 
Lemma 2.2 For any positive integer I, 



2! " 

Proof. This inequality holds for I = 1 obviously. We proceed by induction. 



VW+W _ y/(2/ + 2)! _ y^Q! v /(22 + 2)(2/ + l) 



(2 + 1)! Z!(2 + l) 2! 2 + 1 

< 



VM! yf+f <2 , (22 + 2) = m 
2! 2 + 1 " Z + l 



Our lemma thus follows. jj 
Now we have come to our main result. 

Theorem 2.3 Let h be as in (4) and p be as in Lemma2.1 and the remark following its proof. For 
any positive real number bo, let C = max{-^,8p} and ^0=3^- If f € Ch, m and s is the h spline 
that interpolates f on a subset X of R n , as discussed in the beginning of sectionl, then 



\f(x)-s(x)\ <2 2± *= 1 -ir^ ■ c*c-' v^V3CV5(A')*||/||h (9) 

holds for all x € £1 and < S < So where < A' < 1 is a fixed constant and fl is any subset of 
R n satisfying the condition that for any x £ fl, and any number r with ^ < r < ^(Note that 
gf? < bo), there is a simplex Q, x G Q C Q, with diamQ = r, such that Q contains equally spaced 
centers of degree I, < I < j^g, the centers being contained in X. In particular, if x is fixed, the 
requirement for is only the existence of a simplex Q, x e Q C ft, which satisfies the aforementioned 
properties, and X can be chosen to consist only of the equally spaced centers in Q. 
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The constants (3, a n and Aq are as in (4) and Lemma2.1 and the remark following its proof, 
and the crucial constant X' is given by 

*-(§)* 

and only in some cases mildly depends on the dimension n. The h-norm of f , \\f\\h is defined as in 
the statement following (7). 

Proof. Let p be as in Lemma2.1 and bo be arbitrary positive real number. Fix So = For 
any < 5 < So, it's easily seen that < 3CS < 1. There exists a positive integer I such that 
1 < 31SC < 2. From this we get ^ < I < ^ and < 15 < ^ < b . 

For any x £ Q, let Q be an n-simplex containing x such that Q C £1 and has diameter diamQ = IS. 
Then Theorem4.2 of [8] implies that 

\f(x)-s(x)\< Cl \\f\\ h [ \y-x\ l d\a\(y) (10) 

whenever I > m, where a is any measure supported on X such that 

P(y)dcr(y) = p(x) (11) 



for all polynomials p in PP_ V Here 



whenever I > m. By (8), for all 21 > 2m + 2, 



1 n + ff + l n+1 3-21 



< 2 Ii± ^ ±i 7r :i T i V^^ c ^ c ~'v / ^o(2/3)' (12) 

due to Lemma2.2. 

In order to obtain the bound on \f(x) — s(x)\ as mentioned in the theorem, one only needs to 
find a suitable bound for 

I = ci[ \y-x\ l d\a\(y). 

This can be attained by choosing the measure a in an appropriate way, as the following deduction 
shows. 

Let Y = {xi, ■ ■ ■ ,xn} C X be the set of equally spaced centers of degree I — 1 in Q, where 
TV = dimP[^ 1 . By Lemmal.2, there is a measure a supported on Y such that 

p{y)do-{y) =p(x) 



for all p in P£_ 1 , and 

j q A,m<(T-l 
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Now, 



21 - 3 
I - 1 



2(1- 1) - 1 
Z - 1 



< 



2(1-1) 
I - 1 

1 1 

1 1 



4*-' ^ 6y Stirling's Formula 



■4' 



Therefore, for such <r, we have 
ly-xprf^Ky) < 



n+£+l 2+1 £ _ 

z 4 7T 4 y/na n c 2 c 



< 



< 



< 



< 



n+0+1 .i+l H 

2 4 7r 4 y/na n c 2 c 



n + P + l n + 1 H 

2 4 7r 4 y'nanC 2 c 



n + + 1 ,i + l £ 

2 4 7r 4 ^na n c 2 c 



ii+ff+i n+i . 

2 4 w^* - ,/nanC 2 c 



-/ 



ii+ff+i ,i+i p 

2 3 7r 4 ^Jna n c^c 



^ P ) l (l8) 1 
^(2p) l (lS) 1 



2? - 3 
Z - 1 
1 1 



1 



1 



1 


V^* - 


i) 


1 




V*(f - 


i) 


i 




W*Q - 


i) 



i 



(CIS) where C = max 
i 



2 

3V 



1 < 318C < 2 



2\ ro 
3 



0. Thus \/l — \ ~ whenever <5 is small. Moreover, < I < 



Note that I — > oo as <5 

implies ^ < I < 3C<5 and JfCtf < -77 < V3C^- Thus, for <5 small enough, we get 



hg+1 11+1 



\f(x) - s(x)\ < 2 4 7T 4 ^c^c 'v^o 



1 



40? 



/3CV5(A')*||/IU 



2 ^±|^ ^ VH^cl c" z v^V3CV5(A' )* || /|| h 



2\3C 



where A' = (|) 

Remark. In this theorem although we don't explicitly state that the number 5 represents the 
well-known fill-distance d(Q, Y), it is equivalent to the fill-distance in spirit. This can be seen from 
the inequality < IS < ^ where / denotes the degree of the equally spaced centers Y in Q. Note 
that S — ► iff I — * oo and I — ► oo iff d(Q, Y) — > 0. In this paper we avoid using the term fill-distance 
since the data points are not purely scattered. However the equally spaced centers in a simplex 
are theoretically very easy to handle and in application impose no trouble at all on programming. 
Moreover, in the proof of the theorem, 8 is usually small. Therefore the gap between y/l — 1 and y/l 
can be ignored. This is why we adopt "<" in (9). 

The second thing is about the function space. Although the space Ch, m adopted by us is defined 
by Madych and Nelson in [7] and [5] and is not the same as Wu and Schaback's function space [TU] . 
there is a nice unification theory in llj which is achieved by putting mild restrictions on Wu and 
Schaback's function space. We call it native space. 
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3 Comparison 



The original exponential-type error bound for multiquadric interpolation (See [3] and [5].) is of the 
form 

\f(x)-s(x)\ <C Q \* H/IU, d^O (13) 

, where Co is about the same as the long string of coefficients in (9), A is a constant satisfying 
< A < 1, and d is the fill-distance equivalent to 5 of (9). The definition of A is 



where 



C = max (2//V^e 2 ™ 7 '\ J- 
l 3&o 



c 



with p and c the same as this paper, bo is the side length of a cube, equivalent to the bo of Theorem2.3 
of this paper, and 7„ is defined recurvely by 71 =2, 7„ = 2n(l + 7«-i) if n > 1. The first few j n 
are 

7l = 2, 72 = 12, 73 = 78, 74 = 632, 75 = 6330, • • ■ 

. The fast growth of 7„ as n — > 00 forces e 2 ™ 7 " to grow extremely rapidly, driving A to 1 and making 
the error bound (13) meaningless. 

In our new approach, the crucial constant A' in (9) is defined by 



'I 



— 



where 

f 2 

C = max j— , 8p 

. Note that p is only mildly dependent of dimension n. This can be seen in Lemma2.1 and the 
remark following its proof. A lot of time p, and hence A', is independent of n. For example, when 
(3=1 and n = 1, 2, 3 or 4, A' is completely independent of n. Fortunately, these are exactly the most 
interesting and important cases. For higher dimensions, A' can remain unchanged by increasing p 
in (4) to keep n — (3 < 3. 

Besides this, ^/S in (9) also contributes to the convergence rate of the error bound as 5 — > 0. 
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